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Abstract 

On the basis of angular spectrum representation, a formalism describing paraxial beams propa- 
gating through an isotropic left-handed material (LHM) slab is presented. The treatment allows 
us to introduce the ideas of beam focusing and phase compensation by LHM slab. Because of the 
negative refractive index of LHM slab, the inverse Gouy phase shift and the negative Rayleigh 
length of paraxial Gaussian beam are proposed. It is shown that the phase difference caused by 
the Gouy phase shift in right-handed material (RHM) can be compensated by that caused by the 
inverse Gouy phase shift in LHM. If certain matching conditions are satisfied, the intensity and 
phase distributions at object plane can be completely reconstructed at the image plane. 
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I. INTRODUCTION 



In the late 1960s, Veselago firstly introduced the concept of left-handed material (LHM) 
in which both the permittivity e and the permeability fi are negative Veselago predicted 
that electromagnetic waves incident on a planar interface between a right-handed material 
(RHM) and a LHM will undergo negative refraction. Theoretically, a LHM planar slab can 
act as a lens and focus waves from a point source. Experimentally, the negative refraction 
has been observed by using periodic wires and rings structure 0, [J U, 0, 0| • In the past few 
/ears, negative refractions in photonic crystals QflSQ and anisotropic metamaterials Q 
have also been reported. 



Sars, negative 
,QQ,[l3 



Recently, Pendry extended Veslago's analysis and further predicted that a LHM slab can 
amplify evanescent waves and thus behaves like a perfect lens nj. It is well known that 
in a conventional imaging system the evanescent waves are drastically decayed before they 
reach the image plane. While in a LHM slab system, both the phases of propagating waves 
and the amplitudes of evanescent waves from a near-field object could be restored at its 
image. Therefore, the spatial resolution of the superlens can overcome the diffraction limit 
of conventional imaging systems and reach the subwavelength scale. While great research 



interests were initiated b y t 



raised 



21 



22, 



23, 



24. 



25, 



re revolutionary concept 



26 



27]. 



0,H is, Q, 



hot debates were also 



The main purpose of the present work is to investigate the paraxial beams propagating 
through an isotropic LHM slab. Starting from the representation of plane-wave angular 
spectrum, we derive the propagation of paraxial beams in RHM and LHM. Our formalism 
permits us to introduce ideas for beam focusing and phase compensation of paraxial beams 
by using LHM slab. Because of the negative refractive index, the inverse Gouy phase shift 
and negative Rayleigh length in LHM slab are proposed. As an example, we obtain the 
analytical description for a Gaussian beam propagating through a LHM slab. We find that 
the phase difference caused by the Gouy phase shift in RHM can be compensated by that 
caused by the inverse Gouy phase shift in LHM. 

II. THE PARAXIAL MODEL OF BEAM PROPAGATION 



In this section, we present a brief derivation on paraxial model in RHM and LHM. 
Following the standard procedure, we consider a monochromatic electromagnetic field 
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E(r, t) = i2e[E(r) exp(— iut)} and B(r,t) = _Re[B(r) exp(— iut)] of angular frequency u 
propagating through an isotropic material. The field can be described by Maxwell's equa- 
tions 28( 

^ ^ 9B 

VxE = -i>r- 
VxH = ip 

D = eE, 

B = /iH. (1) 

One can easily find that the wave propagation is only permitted in the medium with e, fi > 
or e, fi < 0. In the former case, E, H and k form a right-handed triplet, while in the latter 
case, E, H and k form a left-handed triplet. The previous Maxwell equations can be 
combined straightforwardly to obtain the well-known equation for the complex amplitude 
of the electric field in RHM or LHM 

V 2 E — V(V • E) + k 2 E = 0, (2) 

where k = n R ^ L u/c, c is the speed of light in vacuum, hr = ^/ErUr and ul = —^/slJIl are 
the refractive index of RHM and LHM, respectively 

Equation (J2J) can be conveniently solved by employing the Fourier transformations, so 
the complex amplitude in RHM and LHM can be conveniently expressed as 

E(rj_, z) = J d 2 k ± E(k ± ) exp[zk ± • r ± + ik z z}. (3) 

Here rj_ = xe x + ye y , kj_ = k x e x + k y e y , and ej is the unit vector in the j-direction. Note 



that k z = 0y ^fj£^o — k 2 ^, a = 1 for RHM and a = — 1 for LHM. The choice of sign ensures 
that power propagates away from the surface to the +z direction. The field ^(kjj In Eq. © 
is related to the boundary distribution of the electric field by means of the relation 

E(kj_) = / dh ± E(r ± , 0) exp[*k ± • r J, (4) 



which is a standard two-dimensional Fourier transform |22|. In fact, after the electric field 
on the plane z = is known, Eq. (jHJ) together with Eq. provides the expression of the 
field in the space z > 0. 

From a mathematical point of view, the approximate paraxial expression for the field can 
be obtained by the expansion of the square root of k z to the first order in |kj_|/A; 0, 13I | - 



which yields 



Efn,*) 



x exp 



ik i z 



E(k ± ). 



(5) 
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Since our attention will be focused on beam propagating along the +z direction, we can 
write 

E(r_i_, z) = A(r_L, z) exp(in R , L k z), (6) 

where the field ^4(r_i_, z) is the slowly varying envelope amplitude which satisfies the parabolic 
equation 







A(r ± ,z) 



0. 



(7) 



dz 2n RtL k 

where Vj_ = d x e x + d y e y . From Eq. (J7J) we can find that the field of paraxial beams in LHM 
can be written in the similar way to that in RHM, while the sign of the refractive index is 
reverse. 



III. THE PROPAGATION OF PARAXIAL GAUSSIAN BEAM 



The previous section outlined the paraxial model for general laser beams propagating in 
RHM and LHM. In this section we shall investigate the analytical description for a beam with 
a boundary Gaussian distribution. This example allows us to describe the new features of 
beam propagation in LHM slab. As shown in Fig.Q the LHM slab in region 2 is surrounded 
by the usual RHM in region 1 and region 3. The beam will pass the interfaces z = a and 
z = a + d before it reaches the image plane z = a + b + d. To be uniform throughout the 
following analysis, we introduce different coordinate transformations z*(i = 1,2,3) in the 
three regions, respectively. 

First we want to explore the field in region 1. Without any loss of generality, we assume 
that the input waist locates at the object plane z — 0. The fundamental Gaussian spectrum 
distribution can be written in the form 



Eifkj 



WQ 

s/2n 



exp 



(8) 



where w is the spot size. By substituting Eq. (JHJ) into Eq. (J3J), the field in the region 1 can 
be written as 




FIG. 1: The mechanisms for paraxial beams propagating through an isotropic LHM slab. The 
LHM slab in region 2 is surrounded by the usual RHM in region 1 and region 3. The solid line and 
the dashed line are the theoretical object and image planes, respectively. 
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(ii) 

(12) 



Here z* = z, z R = n R k wl/2 is the Rayleigh length, Wi{z{) is the beam size and R\(z*) the 
radius of curvature of the wave front. The first term and the second term in Eq. (jl())l are 
the plane wave phase and radial phase, respectively. The third term in Eq. (jl())l denotes the 
Gouy phase is given by $i = — arctan(^/z fi ). 

We are now in a position to calculate the field in region 2. In fact, the field in the first 
boundary can be easily obtained from Eq. Q by choosing z = a. Substituting the field into 
Eq. (JH), the angular spectrum distribution can be obtained as 

k^(n R koWQ + 2ia) 



Eofk 



w 
s/2n 



exp 



(13) 
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For simplicity, we assume that the wave propagates through the boundary without reflection. 
Substituting Eq. (jTHj) into Eq. (JHJ), the field in the LHM slab can be written as 



E 2 (r±,2 2 , 



w 



w 2 {zl) 



exp 



w 2 2 {z* 2 ) 



(14) 
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FIG. 2: Plot of the beam widths in LHM slab with different refractive indices. The solid {q = — 1), 
dash-dot-dotted (ij = —1.2), and dash-dotted (77 = —0.8) lines indicate the beam widths in LHM 
slab, where we assume a = zr and r] = n^jn^. One can find the beams can be focused by LHM 
slab and the focusing waists remain wq. 



Here z\ = z—iX — n^jn^a and Zl = n^kQW^jl is the Rayleigh length in LHM. The beam size 
^2(^2) an d the radius of curvature ^2(^2) are gi ven by Eq. (fTUj) and Eq. (fTTfl. respectively. 
The Gouy phase shift in LHM is given by $2 = — arctan^/z^). 

We note two interesting features of the paraxial field in LHM: First, because of the 
negative index of refraction, the inverse Gouy phase shift and negative Rayleigh length 
should be introduced. Second, the beams can be focused by LHM slab and the focusing 
waists remain wq. For the purpose of illustration, the relevant focusing feature is shown in 
Fig. El We find that the place of focusing waist could be different which depends on the 
choice of the refractive index. 

Finally we want to study the field in region 3. The field in the second boundary can be 
easily obtained from Eq. (|14j) under choosing z = a + d. Substituting the field into Eq. (@J), 



n L k z* 2 + 



n L k r\ 
2R 2 (z*) 



^2 

— arctan — , 

z L 



(15) 
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(17) 



FIG. 3: Numerically computed spatial map of the magnitude of the electric field for Gaussian 
beam propagating through the LHM slab. We assume a = zr, d = 2z R and Ul = — n R . The first 
focusing waist in LHM slab locates at the plane z/zr = 2, while the second focusing waist in region 
3 locates at the plane z/zr = 4. 



the angular spectrum distribution can be written as 



E 3 (k 



W 



± 



exp 



k±(n R riLk WQ + lin^a + 2in R d) 
An R n L k 



Substituting Eq. (fT%|) into Eq. (JHJ), the field in the LHM slab is given by 
E 3 {r ± ,z 3 ) = — — exp 



w 3 {z 3j 



n R k z 3 



n R k r] 



w 3 (z 3 ) = WQ^/lTizlfz^ 



2R 3 (z* 3 ) 



— arctan • 



* 1 
z 3 + 



(18) 

(19) 

(20) 

(21) 
(22) 



Here zl 



o 3 — z — (1 — n R /nL)d. The beam size w 3 (z 3 ) and the radius of curvature R 3 {z^) 
are given by Eq. (J2"T|) and Eq. (J22|) . respectively. The Gouy phase shift is given by $3 = 
— axctait(z 3 / z R ). 

To this end, the fields are determined explicitly in the three regions. Comparing Eq. ()14j) . 
Eq. (fT9*j) with Eq. Q show that the field distributions in region 2 and region 3 still remain 
Gaussian. For the purpose of illustration, the spatial map of the electric fields is plotted in 
Fig. El We can easily find that there exists an internal and an external focus. 



Region 1 



Region 3 




FIG. 4: The phase difference caused by the Gouy phase shift in RHM can be compensated by that 
caused by the inverse Gouy phase shift in LHM slab. The phase fronts of Gaussian beam (solid 
lines) differ from those of a perfect spherical wave (dashed lines). 

IV. BEAM FOCUSING AND PHASE COMPENSATION 



In this section we examine the matching conditions of beam focusing and phase compen- 
sation. First we want to explore the matching condition of focusing. We can easily find the 
place of the focusing waist by choosing z* — 0. We assume the incident beam waist locates 
at plane z = 0. After setting z\ = in Eq. (|16|h we get the first focusing waist in LHM slab 
locates at the plane z = (1 — n^jn^a. Then we substitute z\ = into Eq. (|16J) . we can 
find the second focusing waist in region 3 locates at the plane z — (1 — rifj/n^d. We take 
the image position z = a + <i + 6 to be the place of the second focusing waist. Using this 
criterion, the matching condition for focusing can be written as 



riL^a + b) + tircI = 0. 



(23) 



In addition, the thickness of the LHM slab should satisfy the relation d > —an^jn^ other- 
wise there is neither an internal nor an external focus. 

Next we attempt to investigate the matching condition for phase compensation. It is 
known that an electromagnetic beam propagating through a focus experiences an additional 
7r phase shift with respect to a plane wave. This phase anomaly was discovered by Gouy in 
1890 and has since been referred to as the Gouy phase shift It should be mentioned 
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that there exists an effect of accumulated Gouv phase shift when a beam passing through 



an optical system with positive index 



34, 



35 



361 ] . While in the LHM slab system we expect 



that the phase difference caused by Gouy phase shift can be compensated by that caused 
by the inverse Gouy shift in LHM. In Fig. we plot the distribution of phase fronts in the 
three regions. The phase difference caused by the Gouy phase shift in the three regions are 
given by 

. , a 
A<Pi = — arctan — , 

Zr 

—nia —nib 

A<P 2 = —arctan arctan , 

n R z L n R z L 

b 

A$ 3 = -arctan — . (24) 

zr 

The first and third Equations dictate the phase difference caused by the Gouy shift in region 
1 and region 3, respectively. The second equation denotes the phase difference caused by 
the inverse Gouy phase shift in LHM slab. Subsequent calculations of Eq. (J2"4~j) show 



a b ( —n^a —nib\ 

— arctan arctan h I — arctan arctan = 0. (25) 

zr z r \ n R z L n R z L J 

This means that the phase difference caused by the Gouy phase shift in RHM can be 
compensated by that caused by the inverse Gouy phase shift in LHM slab. Therefore the 
condition for phase compensation can be simply written as 

n R k a + n R k Q b + n L k d = 0. (26) 

The first two terms in Eq. (J26|) are the phase deference caused by the plane wave in RHM. 
The last term in Eq. (f2T)J) is the phase deference caused by the plane wave in LHM slab. 
Finally we discuss the phase difference caused by the radial phase. Following the method 

rj 

outlined by Dumelow et al. [20(, we assume that the beam waist locates at the object plane. 
Then we take the image position to be that where the intensity is a maximum. Evidently, 
from Eqs. ()19|l - (j22|) we can find that the intensity maximum locates at the plane of beam 
waist. The phase difference between the object plane and the image plane is independent 
of radial position, since the phase fronts are flat there. 

The new message in this paper is to prove that the intensity and phase distributions at 
object plane can be completely reconstructed at the image plane by LHM slab. Now an 
interesting question naturally arises: whether the matching conditions of focusing and the 
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phase compensation can be satisfied simultaneously. From Eqs. one can easily find 

that the refractive index, a, b and d should satisfy the matching conditions: 

riL = —tir, a + b = d. (27) 

Under the matching conditions, the reflected waves at the interface between RHM and LHM 
are completely absent. Therefore the intensity and phase distributions at the object plane 
can be completely reconstructed at the image plane. 

Note that the purpose of this paper is to examine beam focusing and phase compensation 
in paraxial regime. The evanescent waves which are claimed to provide the subwavelength 
imaging do not correspond to the problem under study. The paraxial model only deals 
with beams whose transverse dimension is much larger than a wavelength. However, in the 
subwavelength focusing regime a rigorous diffraction theory should be developed. 

V. CONCLUSIONS 

In conclusion, we have investigated the focusing and phase compensation of paraxial 
beams by an isotropic LHM slab. We have introduced the concepts of inverse Gouy phase 
shift and negative Rayleigh length of paraxial beams in LHM. We have shown that the phase 
difference caused by the Gouy phase shift in RHM can be compensated by that caused by 
the inverse Gouy phase shift in LHM slab. If certain matching conditions are satisfied, the 
intensity and phase distributions at object plane can be completely reconstructed at the 
image plane. We expect many potential devices can be constructed based on the paraxial 
beam model discussed above. They can, for example, be used to provide beam focusing, 
phase compensation and image transfer. 
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